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Context. Astronomical observations toward Sagittarius B2(M) as well as other sources with APEX have recently suggested that the 
rest frequency of the / = 1 - transitions of '''CH^ is too low by about 80 MHz. 

Aims. Improved rest frequencies of isotopologs of methylidynium should be derived to support analyses of spectral recording obtained 
with the ongoing Herschel mission or the upcoming SOFIA. 

Methods. Laboratory electronic spectra of four isotopologs of CH+ have been subjected to one global least-squares fit. Laboratory 
data for the 7=1-0 ground state rotational transitions of CH^, "CH^, and CD^, which became available during the refereeing 
process, have been included in the fit as well. 

Results. An accurate set of spectroscopic parameters has been obtained together with equilibrium bond lengths and accurate rest 
frequencies for six CH+ isotopologs: CH+, "CH+, "CD^, CD+, ''*CH+, and CT+. 

Conclusions. The present data will be useful for the analyses of Herschel or SOFIA observations of methylidynium isotopic species. 
Key words, molecular data - radio lines: ISM - Submillimeter: ISM - ISM: molecules 
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1. Introduction 

The first three molecules detected in space in the 1930s were 
methylidynium, CH+, methylidyne, CH, and cyanogen, CN. 
iDouglas & Herzberg (1941) identified three absorption lines, 
which had been observed toward several stars, as belonging to 
the A - X 'S^ electronic transition of CH^. They were as- 
signed to the R(0) {= J - 1-0) transitions of the y = 0-0, 1 - 0, 
and 2-0 vibrational bands. 

The 7 = 2 - 1, 3 - 2, and 4-3 r otational transitions 
were identified by ICernicharo et al.l (Il997h in spectra recorded 
with the Long Wavelength Spectrometer of the Infrared Space 
Observatory toward the planetary nebula NGC 7027. The obser- 
vation of the J =1-0 rotational transitions at 835137.5 MHz 
(lAmanoll2010|) is hampered by a telluric line of molecular oxy- 
gen close bv (iDrouin et al.ll2010l) and thus requires observations 
from space, e.g. with the recently launched Herschel satellite 
or with the Stratospheric Observatory For Infrared Astronomy 
(SOFIA). 

Recently, iFalgarone et al.l (l2005l) reported on the obser- 
vation of the J - 1-0 transition of '^CH^ in absorp- 
tion toward the star-forming region GlO.6-0.4. They deduced 
830132 (3) MHz as the prefered rest frequency based on scal- 
ing from the CH^ laboratory rest frequency. Alternative rest fre- 
quencies of 830107 (1) or 830193 (4) MHz were excluded. Very 
recently, Menten et al. (2010) carried out an absorption study 
with the Atacama Pathfinder Experiment (APEX) of light hy- 
dride species toward the evolved massive star-forming region 
Sagittarius (Sgr for short) B2(M). They detected low-lying rota- 
tional transitions of '^CH+, H^^Cl, H^'^Cl, and, for the first time, 
of SH^. Similar studies have been carried out toward additional 
sources (F. Wyrowski, private communication). A comparison of 
the line profiles among these species together with observations 



of additional species at lower frequencies suggested a rest fre- 
quency of about 830210 MHz, approx imately 80 MHz highe r 
than the prefered rest frequency from IFalgarone et al.l ( l2005h . 
but in reasonable agreement with an alternative rest frequency 
of 830193 (4) MHz they had dismissed. 



iPearson & DrouinI (l2006l) obtained a value of 
835078.950 (75) MHz as the 7 = 1-0 transition fre- 
quency of CH^, about 5 MHz lower than they deduced from 
analyses of the A '11 - X '2^ electronic spectra of four 
i sotopologs (Carrin gton & Ramsay 1982j .Bembenek ,199731 : 
IBembenek et al.|[l987l:rBembeneklll997bl) . 

During the refereeing process of this manuscript I received 
7=1-0 transition frequ encies forCH+, '^CH+, and CD+, 
very recently obtained by Amand (|2010|) . While the latter two 
frequencies (~830216 and 453522 MHz) were rather close to 
predictions made from the spectroscopic para meters of the re - 
spective isotopolog ( Bembenek 1997a; Bem benek et alj[T987l) . 
the CH+ frequency (835137.5 MHz) was almost 40 MHz higher 
than that deduced from a more recent analysis of the CH^ 
electronic spectrum (Hakalla et al. 2006), but almost 70 MHz 
higher than measurements by Pearson & DrouinI ([2006) indi- 
cated. Because no further rotational data and no infrared tran- 
sitions have been reported thus far, determinations of methyli- 
dynium rest frequencies still rely quite heavily on the electronic 
spectra. 



In the current work, I present a combined reanalysis of the 
A 'n - X '2^ electronic spectra of four isotopic species of CH^ 
together with the available rotational transitions to derive reli- 
able rest frequencies for CH^ isotopologs. 
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2. Analysis and discussion 

The ro- vibrational energy levels of a diatomic molecule AB can 
be represented by the Dunham expression 



E{v, j)ih = Yj ^'M + layj^ij + ly 



(1) 



where the Yij are the Dunham parameters. In electronic states 
different from S states, i.e. in states with orbital angular momen- 
tum A > 0, the expansion in J(J + 1) is commonly replaced by 
an expansion in J(J+ 1)-A^, see e.g. (IWatsonll980h . Obviously, 
the two are the same for S states with A = 0, such as the ground 
electronic state of CH^. The A state is a '11 electronic state with 
A = 1, and the expansion has been carried out in J(J + 1) - 1; an 
expansion in J( J +1) can be found also, see e.g. the case of the 

radical BrO (iDrouin etalJl200lh . 

The energy levels in a 'll electronic state are modified by a 
/l-doubling which can be expressed as 



E'\v, J)/h = +J(J + l)/2 ^ qijiv + 1 /2yjj(J + ly. 



(2) 



-2-(/+2;)/2 



where the /l-do ubUng parameters qij scale with fi^ 
(lBrownetalJll979|). 

Watsonl(ll973Lll98d) has shown that several isotopic species 



of AB can be fit jointly by constraining the Yjj to 



nipA 



1 + 



Ma 



m A^^ 



Mr 



-(/+2;)/2 



(3) 



where Uij is isotope invariant, m^, is the mass of the elec- 
tron, He - Ma Mb /(Ma + Mb - nwg) is the charge-corrected 
reduced mass of AB and n is the charge of AB, Ma is the 
mass of atom A, and Afj is a Born-Oppenheimer breakdown 

(BOB) term. Uijii'^'^^^^^^AfjinelMA is sometimes abbreviated as 
6fj. Moreover, it is noteworthy that both Af. and 6fj are defined 
negatively in some papers. Obviously, zl^ and 5? are defined 
equivalently. 

A specific transition of the A '11 - X 'S^ electronic spec- 
trum of CH^ is well described by the isotopolog, the vibrational 
quantum numbers in the upper and lower electronic states, the 
lower state rotational quantum number J and whether the tran- 
sition belongs to the R, Q, or f -branch in which the upper state 
rotational quantum number is 7 -i- 1, 7, and 7-1, respectively; J 
is commonly given in parentheses. 

None of the investigations of the A - X electronic spectrum 
of CH^ and its isotopologs provided uncertainties for individual 
lines, only estimates of average precision and accuracy. In the 
present work the initial uncertainties were estimated for each 
isotopolog and each vibronic band separately based on aver- 
age residuals between observed frequencies and those calculated 
from the spectroscopic parameters. During the analysis, the un- 
certainties were adjusted to yield an rms error very close to the 
theoretical value of 1.0, usually slightly lower. Transitions with 
residuals between observed and calculated frequencies much 
higher than the average were omitted from the fit. These may 
have been weak lines with very low or very high values of J 
or accidentally overlapped lines. For each isotopolog and each 
vibronic band the finally ascribed uncertainty will be given be- 
low in parentheses in units of 0.001 cm ' together with the lines 
from the original literature, which have been omitted from the 
fit. Lines already discarded in the original analyses were usually 



Table 1. Experimentally determined spectroscopic parameters" 
(cm^', MHz) of methyhdynium, '2CH+ in its X 'E+ and A 'n 
electronic states from a global fit involving four isotopologs. 



Parameter 



Value 





2 860.387 (22) 




-0.147 (18) 




-2.672 (13) 




-59.369 7 (61) 




0.05 1 (62) 




0.225 46 (106) 




427 360.68 (121) 


UniU^'A^.m, /Mr 


-155.80 (21) 


^ u I " (} 1 ^ / n 


-2 146.22 (177) 




-14917.7 (57) 


Unl^'^'^A^itnJMH 


58.6 (27) 


Y21 


91.8(41) 


Y31 


-3.84 (66) 




-42.468 (21) 


UoiiJ-^A^^mJMH 


0.563 (20) 


Y12 


0.960 (34) 


Y22 


-0.045 (10) 


Yo3 X 103 


4.264 (86) 


Kb X 10^ 


-0.726 (53) 


Yo4 X 10*'' 


-0.4 


cC^c) 


1.091 (28) 


A 


'n 




23 639.326 (23) 


UooA'^^mJMcillf 


-19.620 (24) 




23 786.537 (24) 


UooA'iamJMciUt 


-26.436 (26) 


t/ioyu;"'(H) 


1 864.379 (14) 


Uiatf,"^A%mJMc(n) 


0.106(12) 




1 367.473 (17) 




-0.180(19) 


YioQl) 


-115.8847 (48) 


Y20(O) 


-60.815 4(47) 


Y3oQi) 


2.639 48 (85) 


^30(0) 


0.776 92 (78) 




356717.5 (63) 




-9.2 (53) 




-316.4 (24) 




-27 562.7 (70) 


Unf^c^'^A^^mJMn 


24.06 (147) 


Y21 


-634.8 (49) 


^31 


137.76 (95) 


Ym. 


-58.166 (26) 


Yn 


-2.502 (61) 


Y22 


-0.252 (44) 


Y32 


0.149 3 (87) 


Yo3 X 103 


-0.876 (84) 


713 X 10^ 


-0.911 (63) 


qoo 


1234.25 (108) 


910 


-96.53 (82) 


901 


-0.643 6 (98) 


911 


-0.020 3 (79) 



" The vibrational parameters Yjo as well as related parameters are given 
in cm"'; all other Yjj as well as related parameters are given in MHz. 
Numbers in parentheses are one standard deviation in units of the least 
significant figures. 

* Determined with respect to the lowest allowed energy levels of each 
electronic state. 
Estimated, see text. 
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omitted here as well; these lines will not be mentioned specifi- 
cally. 

The fit was s ta rted with t he CH^ da ta from 
ICarrington & Ramsavl (Il982h and iHakalla et all (l2006h . 
and the initial parameters were taken from the latter work. As 
the Yoi got s maller by facto r s of ~ -10000, Yq4 was estimated 
as -0.4 Hz. IHakalla et all ( 12006 ') measured the y = 0-0 
(5 (xO.OOl cm"'), j;(17)), - 1 (5 -) and 2 - 1 (9, P(14)) 
bands, while ICarrington & Ramsa \l (ll982h obtained data for the 
i; = - (5, P(6)), - 1 (6, 2(14)), 1 - (12, R(3)), 1 - 1 
(7, P(7,8)), 1 - 2 (6, P(7)), 1 - 3 (4, R(5)), 2 - 1 (5, -), and 
3 - 1 (7, Q(9), P(2)) bands. The CH+ 7=1-0 rotational 
transition predicted from this data set is 8350 93.42 + 5.25 MHz, 
almost exactly 10 MHz higher than what iPearson & DrouirJ 
daooe") obtained and about 14.5 MHz higher than the line these 
authors measured in the laboratory. Inclusion of data from 
three additional isotopologs, as described below, shifted the 
prediction to somewhat higher frequencies, ~835112 MHz 
being the highest as long the the weights of the CH^ data were 
not reduced. 

Scaling these parameters directly with fj.c yields 
830146.97 MHz for the ' ^CH+, 7=1- transition. In 
contrast, predictions from lBembenei3 (Il997ah yield a value 
approximately 80 MHz higher, suggesting a large BOB term 
for Foi will be necessary at least in the ground electronic state. 
Inclusion of data from the i> = - (4, -) band confirmed this. 
These data also required a BOB term for Fqi as well as for the 
term value T of the excited electronic A state, though the latter 
value appeared to be too high to be considered a BOB term in 
the usual sense. The 7=1-0 ground state rotational transition 
was now predicted at 830206.82 ± 5.85 MHz, 60 MHz higher 
than from direct scaling as a result of large effects caused by 
the breakdown of the Born-Oppenheimer approximation. The 
y = - 1 (4, -), 2 - (8, Ri^)) and 2 - 1 (12, -) bands were 
included without much difficulty. In the v - 1 - (6, -) band, 
in which already several lines were not used in the original fits, 
the ^(5,9) and f (3,4, 10) lines were omitted in addition in the 
present work. In the v = 1 - 1 (9, -) band, R(6,l) are in fact 
7?(7,8), similarly the 7 values have to be revised one up for 
g(7,8). Moreover, R(1,9,10) were omitted from the fit. With 
more vibronic bands included, a correction term for Fio was 
also employed in the fit. 

The CD+ data from iBembenek et all (1198 7') proved to be 
much more difficult to reproduce. Not only were BOB terms re- 
quired for the Yoi and Fio, as for '^CH^, but also for Y20, Fi 1 , and 
Y02 in both electronic states. The correction term to T turned out 
to be very large again. The y = - (4, -), - 1 (4, ^(12)), 1 - 
(3, Rill)), P(2,6,12,13)), 1 - 2 (9, -), 2 - 1 (6, -), and 3 - 1 (5, 
Q{1), P(6,ll)) bands were used in the fit. The v - I - 3 ban d 
was measured but not used in the fit bv lBembenek et all (Il987h . 
it is on average about 0.09 cm"' too low. In addition, y = 2 - 
was also not used in the present fit as it was about 0.03 cm"' too 
low. Its inclusion in the fit would have afforded both y = 2 - 1 
and y = 3 - 1 to be omitted f rom the fit. 

Lastly, '^CD^ data from iBembenekl (Il997bl) were included 
in the fit, namely the t; = - (6, f (12,13,14)) and 1 - (4, -) 
bands. The band origins were not quite compatible with the other 
data. In combination with the considerable corrections required 
for T, this was interpreted as a sign of the perturbation of the 
A electronic state. Therefore, separate term values as well as y,o 
were used for CH^ and CD^. In addition, different carbon BOB 
terms were required for T and Fiq. 

The referee of this manuscript pointed out new laboratory 
measurements of the 7=1-0 rotational transitions of CH^, 



Table 2. Derived spectroscopic parameters" (cm"', MHz, unit- 
less) of methylidynium, '^CH^ in its X 'S^ and A 'll electronic 
states from a global fit involving four isotopologs. 



Parameter 



Value 



^10 



A^ 



A'^ 
A^ 

5^11 
A^ 

Ykl 
A^ 



2 857, 
-1, 
-1, 
-59, 
-1, 
425 058, 
-7, 
-9, 
-14 859, 
-7, 
-41, 
-24, 



569 2 (92) 
120(137) 
715 9 (86) 
318 7 (60) 
578 (188) 
66 (248) 
9749 (105) 
226 3 (76) 
0(67) 
222 (33) 
905 (20) 
36 (86) 



A 'n 



r(H)'' 

r(D)'' 

^^(H) 

^m(H) 
i'io(D) 

^m(D) 

A'^ 

A^ 
^01 



^11 



23 619, 
23 760, 
1 864, 
1, 

1367, 
-2, 

356391, 
-0, 
-1, 

-27 538, 
-1, 



707 (5) 
102 (5) 
485 (8) 
244 (145) 
294 (8) 
876 (306) 
8(29) 
56 (33) 
630 (12) 
7(71) 
604 (98) 



" The electronic term values T , and the vibrational parameters F/o are 
given in cm"'; all other Yij are given in MHz; the zl^ are unitless. 
Numbers in parentheses are one standard deviation in units of the least 
significant figures. 

* Determined with respect to the lowest allowed energy levels of each 
electronic state. 



'^C H"^, and CD"^ w hich I received soon thereafter from the au- 
thor (lAmanol2010h . The '^CH^ and CD^ fr equencies were clos e 
to the prediction s based on parameters from lBembenekl(ll997al) ; 
Bemben ek et al ] (119871) as well as predictions from present fits 
which only included data from the electronic spectra. In con- 
trast, the CH^ frequency was higher than predictions from the 
cuii-ent fits by between 25 and 45 MHz. The CH+ and '^CH+ 
lines showed large Zee man splitting in a magnetic field, which, 
as already mentioned by Amano (2010), requires a high g y-value 
for a 'Z^ molecule. This in turn leads to high Aoi BOB terms for 
CH^ in the ground electronic state, see below. The large gy-value 
is also compatible with the large '^C nuclear spin-rotation split- 
ting of 1.6 M Hz in the 7=1 - transition of '^CH+, as already 
pointed out bv lAman ^ (l20T0l) . 

The fit easily accomodated these rotational data; the rms er- 
ror of the fit only slightly increased from 0.944 to 0.991. This 
rather small deterioration of the fit may indicate unfavorable cor- 
relation of the parameters, which is reduced by the inclusion of 
the pure rotational data in the fit. Weighting two v ibronic bands 
of CH + each from ICarrington & Ramsavl dl982l) : iH^alla et al.1 
(l2006h slightly lower (< 25 %) yielded rms errors of all vibronic 
bands very close to 1.0 and an overall rms error of 0.945. The 
inclusion of the rotational data caused the '^C BOB correction 
value to f/oiA*"' in the ground electronic state to be reduced by 
almost half, similar changes in magnitude occured in the excited 
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Table 3. Measured (lAmandbolOh and calculated low- J rota- 
tional transition frequencies" (MHz) of six isotopic species of 
CH+. 



Measured 



Calculated 





7=1-0 


7=1-0 


7 = 2-1 


CH+ 


835 137.504 (20) 


835 137.503 (20) 


1669 281.3 (3) 


"CH+ 


830 216.095 (30)* 


830216.096 (22) 


1659450.3 (3) 


CD+ 


453 521.851 (20) 


453 521.851 (20) 


906 752.0(1) 






448 539.258 (26) 


896 793.2(1) 






826012.748 (45) 


1651053.6 (3) 


CT+ 




325 443.364 (70) 


650736.8 (1) 



" Numbers in parentheses are one standard deviation in units of the least 
significant figures; see the text for discussion of uncertainties. 
' Value corrected for hyperfine splitting; two lines were measured, see 
text, each accurate to 30 kHz. 



A state '^C BOB term. These changes were mostly compensated 
by changes in UmfT^ , to a much lesser extent by the BOB term 
for H. Changes considerably outside the uncertainties also oc- 
cured for t/o2jU^ and Fo2 in the ^ and A state, respectively. All 
other changes were mostly well within the uncertainties. 

The final spectroscopic parameters are given in Table [T] fur- 
ther spectroscopic parameters, which were derived from these 
parameters, are given in Table |2l Predictions of the 7=1-0 
and 2-1 transitions of six isotopologs are given in Table |3] to- 
gether with the measured ones; these should be useful for astro- 
nomical observations as well as for laboratory spectroscopic in- 
vestigations. The three laboratory rotational transitions ( Amano 
12010 ') determine three parameters related to Uq\. In addition, 
the weight of each line in the fit is much larger than that of a 
rovibronic transition. Therefore, it is not surprising that in these 
cases predicted frequencies and their uncertainties are essen- 
tially identical to laboratory values (Amano 2010). More ex- 
tensive predictions, including those of rovibrational spectra, as 
well as details of the fit file will be available in the Cologne 
Database fo r Molecular Spectroscopy, CDMSQ (Mull er et aP 
I2OOILI2OO5I) . Extrapolations beyond 7=1-0 should be viewed 
with some caution, meaning the actually measured transition fre- 
quency may be more than three times the predicted uncertainty 
away from the calculated position, especially for the CH^ main 
isotopolog because the rotational transition was rather far away 
from prediction absed on electronic spectra. Predictions with 
7 > 10 should be vie wed with great caution . 

In cont rast to e.g. Carrington & Ramsavl (Il982h . but in agree- 
ment with iHakalla et all (12006'). the excited A electronic state 
could be reproduced in a rather satisfactorily manner using 
power series in y + 1/2 besides 7(7 + 1) - 1 for each isotopolog 
individually and, with the exception of the term value and the 
vibrational spacings, also for the combined data set of four iso- 
topic species. The fact that separate parameters had to be used 
for CH^ and CD^ may be explained by a perturbation of the A 
state. Essentially all parameters Yij with y > or related param- 



eters Uij^, 



-('■+2»/2 



in Table[T]as well as the BOB terms in Table|2] 



appear to be reasonable. 

The Yij or related parameters for the X ground electronic 
state in Table[T|again appear rather reasonable. The BOB terms 
AjQ in Table |2l are negative and their ma g nitudes are normal 
(Watson' 119731 119801: iDrouin et al.1 1200 U iMiiller et al.l 120051 
12007) . The BOB terms however, are rather large, but of 



website: |http://www.astro.uni-koeln.de/cdms/| 



similar magnitude. IWatsonI (Il973h has shown that A{)\ has one 
contribution gjiiie/mp, where gj is the rotational g-factor and 
nip is the mass of the proton, which may get big for low-lying 
electronic states of the same multiplicity, i.e. singlet states in 
the case of CH^. Because Yn is the first vibrational correction 
term to Fqi it is probably reasonable that zl^j is large in mag- 
nitude also. iGauss & Puzzarinil (|2010') have shown that Aqi can 
be calculated ab initio and presented resul ts which were com- 
pared with experimental data for SiS (Mu ller et al.ll2007h . HCl, 
and HE. Calculations for CH+ will be very welcome. 

Equilibrium structural parameters of 112.77543 (16) and 
123.45370 (106) pm have been derived for the X and A electronic 
state, respectively, in the BO approximation which are very dif- 
ferent from 1 13.08026 (33) and 123.50991 (50) pm obtained for 
the CH^ isotopolog under consideration of the breakdown of 
the BO approximation. The latter values agree fa irly well with 
1 13.08843 (30) and 123.5053 (37) pm obtained bv lHakalla et al.l 
(l2006») . 



3. Conclusion 

The A 'n - X 'S^ electronic spectra of four isotopologs plus the 
rotational data for three species have been reproduced rather well 
with one set of spectroscopic parameters and with only compar- 
atively few transitions omitted from the fit. Considering that the 
A electronic state has been judged to be perturbed by an uniden- 
tified state, the spectroscopic parameters show few peculiarities. 

The predicted 7=1-0 transition frequencies of the three 
remaining isotopologs as well as higher-7 predictions for all 
species will be useful for future astronomical observations or 
laboratory spectroscopic investigations. The data for '"^CH^ and 
CT+ may be useful for determining if ''*C or T play a significant 
role in certain (circumstellar) environments. 
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